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SURFACE roughness EFFECT ON FINITE OR JOURNAL BEARINGS 
B. C. Mfijumdar* and B. J. Haimock 
National Aeronautics and Spate Auminlslrat ij)ii 
Lewis Research Center 
Cleveland, Ohio 4413S 
SUMMARY 

A theoretical study of the performance of finite oil journal bearings 
is made considering the surface roughness effect. The total load-supporting 
ability under such a condition derives from the hydrodynamic as well as 
asperity contact pressure. These two components of load are calculated 
separately. The average Reynolds equation for partially lubricated surfaces 
is used to evaluate hydrodynamic pressure. An analytical expression for 
average film thickness is obtained and introduced to modify the average 
Reynolds equation. The resulting differential equation is then solved 
numerically by finite difference methods for mean hydrodynamic pressure, 
which in turn gives the hydrodynamic load. Assuming the surface height 
distribution as Gaussian, the asperity contact pressure is found. The 
effect of surface roughness parameter, surface pattern, eccentricity ratio, 
and length-to-di ameter ratio on hydrodynamic load and on side leakage is 
investigated. It is shown that hydrodynamic load increases with increasing 
surface roughness when both journal and bearing surfaces have identical 
roughness structure or when the journal only has a rough surface. The trend 
of hydrodynamic load is reversed if the journal surface is smooth ana the 
bearing surface is rough. 


‘National Research Council - NASA Lewis Research Center Associate. 



INTRODUCTION 


The first step in gaining insight into the lubrication of solid sur- 
faces is to examine the surface profile or topography. Smooth surfaces are 
not flat on an atomic scale. Tlie roughnesses of manufactured surfaces used 
in lubrication arc between IxlO” and 20x10”° m, whereas typical atomic 
diameters are between 1x10“^^ and 10x10”^^ m. Even a highly polished 
surface, when examined microscopically or with a prof i lometer, has an ir- 
regular nature. The surface consists of high and low spots. The high 
spots, or protuberances, are also called asperities. 

When the bearing surfaces are not perfectly smooth, there is a possi- 
bility of asperity contact. The contact pressure thus developed on asperi- 
ties will carry a part of the applied load. As the surfaces will have 
relative motion under lubrication, hydrodynamic pressure will also be gen- 
erated as a result of wedge action and the viscosity of the lubricant. 

A few theories for predicting the surface contact pressure of nominally flat 
surfaces are available [1-4]. There are two main approaches [5,6] for cal- 
culating the hydrodynamic load of partially lubricated surfaces. Patir and 
Cheng [5] used the flow simulation method of a randomly generated rough sur- 
face with known statistical properties over the surface area. T0nder [6] 
studied the lubrication of a rough surface by a Monte Carlo method. In an 
earlier report the authors [7] have estimated the mean hydrodynamic and 
contact loads of rough surface plane sliders and infinitely long journal 
bearings with both surfaces having identical roughness structures. As both 
the surfaces had the same roughness structure, the shear flow effect was 
neglected in [7]. 

The objective of the present paper is to study the surface roughness 
effect on finite journal bearings where the two surfaces have different 
roughness structures. The present work uses a flow model similar to that of 



Patir and Cheng [b], but the average gap height is obtained analytically. 
This is then introduced to modify the governing differential equation. The 
contact load is computed by using a model given by Greenwood and Tripp [4j. 
The partial differential equation for mean hydrodynamic pressure is solved 
numerically by finite difference methods satisfying Reynolds boundary condi- 
tions. The hydrodynamic load, attitude angle, and side leakage are found 
fn various length-to-diameter ratios, eccentricity ratios, roughness 
parameters, surface patterns, and variance ratios. The effect of these 
design parameters on hydrodynamic load and side leakage is discussed. 

NOTATION 

C nominal radial clearance [m] 

D journal diameter [m] 

E modulus of elasticity of bearing material [N/m^] 

E' composite modulus of elasticity of two surfaces [N/m^], 



e 


h,H 

K 


L 


n 


P.P 


q.Q 

R 


eccentricity [m] 

nominal film thickness [m], H = h/C 

average film thickness [m] 

a constant 

bearing length [m] 

number of asperities per unit area 

mean hydrodynamic pressure [N/m^], P « pC^/6nuR 

p 

contact pressure [N/m ] 

side leakage [m^/s], Q - (L/D)q/CRu 

journal radius [m] 
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f 


velocity of surface [m/s], u ■ Uj, • 0 

•j 2 

variance ratio, V^j - (oj/o) , ■ {o^lo} 

contact load [N], - Wj./LDL* 

V 2 

hydrodynamic load [N], Wj^ * w^C /6nLR u 
coordinates [m], Y « y/(L/2) 
mean radius of curvature of asperities [mj 
surface pattern parameter 
cocnbined roughness [m], 6 ■ 6. + 
roughness ^rp'itude of journal and bearing surfaces 
eccentricity ratio, c - e/C 
absolute viscosity of lubricant [N s/m^] 
angular coordinate [rad], e « x/R, • angular coordinate 
where film breal.s 
roughness parameter, A » C/o 
Poisson's ratio 

standard deviation of combined roughness, o ■ 
standard deviation of journal and bearing surfaces 
shear flow factor associated with two surfaces 
pressure flow factors 

shear flow factor associated with single surface 
attitude angle [rad] 

bearing surface 
journal surface 

THEORETICAL ANALYSIS 


The average Reynolds equation for partial hydrodynamic lubrication has 
been derived by Patir and Cheng [5], and it can be written after neglecting 
the local squeeze term for a journal bearing (fig. 1) as 
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( 1 ) 


fr (‘x fe? ^ (*y TO ^) ■ "V 


ax 




•*! 

ax 


where hj is given by 


ij - ^ (h ♦ «)f(«)d« 


( 2 ) 


and f(6) is the probability density function of combined roughness a. 

The flow factors t and t> will approach 1 as h/o approaches 

X / 

*, whereas will be equal to 0 for a large value of h/g. The flow 

factors b , P , and b not only depend on h/o, but are also func- 
X y s 

tions of the statistical properties of the frequency density of roughness 

heights and the directional properties of asperities. The flow factors 

b , b , and b. were obtained by Patir and Cheng [b] through flow 
X y s 

s’l.iulation of a rough surface having Gaussian distribution of surface 
height. These factors are used in the present calculation. The average 
gan height h^ is calculated iti the followincj way: 

for a Gaussian distribution, the normal probability density function 
uf 6 is 

f(«) = —i— (3) 


Substituting equation (3) into equation (2) gives 


After performing integration. 


i- 


«)e-‘ '2° d6 


(«) 




2 2 

Wile 
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Differentiating hj with respect to x, we get 



The flow factors and Py are given by 

Pjj - 1 - for t < 1 

and 

*y (^ ’') • *« (^ t ) 

vihere c^ and r are constants which can be found in [8j and y is de- 
fineo as the ratio of lengths at which autocorrelation functions of the x 
and y profiles reduce to 50 percent of the initial value. This y can be 
tiunight ot as the lertgth-to-widtii rutU) of a representative asperity. Per 
the definition, purely transverse, isotropic, and longitudinal roughness 
patterns correspond to y of 0, 1, and ®, respectively. 


Knowing the flow factors and fcj. and assuming u^ ■ 0 and Uj - u, 
we can write equation (1) for a constant n as 



When h/o approaches a large value (i.e., for smooth surfaces), equation 
(6) reduces to the classical two-dimensional Reynolds equation. It has been 
found that if h/o ^ 6, the surface roughness effect can be neglected. 

Hence a surface having h/o * 6 corresponds to a smooth surface. 
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Patir and Cheng [B] have shown that the shear flow factor 

represented by the following relationship, which is valid when 

^ - (V . - V . )» 

' rj rb' s 


can be 


where 


and 



Substituting 



1 - V^j in the above expression of 
♦s- - 


»5 - A 


.( 5 ) 




for ^ < S 
a — 


^ -0.25(h/o) h ^ ^ 


where Aj^, A^, oj, a^, and are constants that can be found in 

[ 8 ]. 

Depending on the roughness configuration, can be positive, nega- 
tive, or zero. To understand the physical significance of the sign of 
let us assume surface j is moving and surface b is stationary (fig. 2). 
The additional flow transport due to the combined effect of roughness and 
sliding is (u./2)o^ . If surface j is rough, the fluid carried in the 

J ^ 

valleys results in an additional flow transport, thereby meaning a positive 
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On the other hand, If surface j is smooth, the asperities of the 
stationary rough surface act as barriers 1n restricting the flow. This 
gives a decreased flow as revealed by a negative 6^. 

In the present case, when V _ > 1, ^ Is positive, which means a 

• J W 

rough moving surface and a smooth stationary surface. In other words the 
journal has a rough surface and the bearing a smooth one. When V . ■ 0, 

* J 

the journal will have a smooth surface and the bearing a rough one. 

When ■ 0.5, Oj ■ This means both the journal and bearing 
surfaces have equal roughness. In this particular case, there Is no 
dj effect. 

For the journal bearing shown in figure 1, when the journal is rotatit-.g 
with surtace velocity u, it carries a hydrodynamic load w^ in addition 
to asperity contact load. Assuming h and are functions of x only 
(implying no misalignment) and using the following substitutions 



Dimensionless mean film thickness H is given by 


H » 1 + f cos e 

Thus aH/ae « >c sin e, and equation (7) reduces to 


H 





/P 


3 >P *^x 


•j[‘-"'fe)]‘- 


‘ »'" *1 ■ ’> 


The boundary conditions are 


p . 0 at e - 0 


P • — *0 at 9 ■ 9, 

39 £ 


>0 at I ■ 0 


P - <0 at Y - 1 


where is the anoular coorainate at which the film cavitates. These 

conditions are known as Reynolds boundary conditions. 

For a particular set of values of L/D, e, y (i.e., b , ^ , ana 

X y 

4 ,). V and A, equation (8) )s solved numerically by finite differ- 

^ V 

ence methods with successive overrelaxation factor satisfying the above 
boundary conditions. The numerical calculation is done with an accuracy of 
O.Ol percent of the difference of integrated pressures of two successive 


iterations. 


Calculation of hydrodynamic Load 


Knowing the hydrodynamic pressure distribution, we calculate the two 


components of hydrodynamic load 


r 112 /• 9. 

/ I ‘ 

-ro Jq 


pR cos 9 d9 dy 
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or In dimensionless form 


W 


X 


and 




pfi sin • de d> 






cos • de dY 



> 



sin 6 de dY 


J 


where 


(11) 






and 


“v - 

^ 6nuLR^ 

The integrations of equation (11) are performed numerically by Simpson's 
rule. The total hydrodynamic load Wf, Is 


\ 




7 ? 

y 


( 12 ) 


where 


u 

•^h 



6nuR^ 


The attitude angle ♦ is 
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‘■“"■'ft) 

Calculation of Side Leakage 
The bearing side leakage can be calculated from 

^ ly 

•'o 

A >P 



(13J 


(14) 


where 


Q 


(fe)- 

‘TRir 


(15) 


Iho differentiation aP/aY of equation (14) is obtained numerically by the 
3-point backward difference rule, and then the integration is performed 
numerically by Simpson’s rule. 

Calculati on of Asperity Contact Load 
By using a Ga issian distribution of asperity height, the contact load 
can be evaluated from nominal contact pressure. When both surfaces are 
rough, the nominal contact pressure given by Greenwood and Tripp [4] is 


P 


c 


KC’F 


5/2 


(lb) 


where 


K 


«(nao} 
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and 


V77 4 

The function calculated In [4j for various A ana tabu- 

lated, and K can be calculated for a particular o/a from the above 
expression. 

The contact load ii 

"c ■ '■“'’c 


or 


W 


c 



( 17 ) 


where 


w 


c 


njT' 


U«) 


The two dimensionless loads and are definet in such a way 
that they cannot be simply added. In the following section a numerical 
example is taken to show the procedure for calculating these loads. 

Example 

An oil journal bearing is operating under the following conditions: 
Diameter of jo-rnal 0 ■ 50 mm 
Length of bearing L • 50 nri 
Radial clearance C ■ 0.01^5 nm 
Eccentricity ratio c ■ O.ii 
Journal speed u - 5 m/s 
Absolute viscosity oil « ■ 0.01 N s/m^ 
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Assume E* ■ 2.2x10^^ N/m^ and o/B - 0.01. Find the contact load and 
the hydrodynamic load. 

Although will be different for different h, F ^^2 


5/2 

h . h^.,^ is used. 


^in - 


Substituting the values of C and c, 


hijiii «= 0.0050 min 


The contact load is 


Substituting the values of 


0 



• L0KE-F5,2 
L and taking 


K 


0.003 give 


« 165x10^ N 

Assume a roughness parameter A of 1. The function corresponding 

to A ■ 1 is 0.08056 [4]. Thus the contact load w^ is 133 kN. 

It is further assumed that both journal and bearing surfaces have the 
same roughness structure and that the surface roughness pattern is iso- 
tropic. From figure 3 the dimensionless hydrodynamic load for the above 
bearing is « 2.6. 

Hence the hydrodynamic load is given by 


6nLR u W. 


Substituting the above data in this expression yields 

» 39 kN 

The ratio of contact to hydrodynamic load is 


w 

c 



3.41 
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RESULTS AND DISCUSSION 

The hydrodynamic load, attitude angle, and side leakage in dimension- 
less form are computed for various length-to-diameter ratios (0.5, 1.0, and 
?.0), eccentricity ratios (0.2, 0.4, 0.6, and 0.8), roughness parameters 
(1, 2, 3, 4, and 6), surface pattern parameters (1/6, 1, and G), and vari- 
ance ratios (0, 0.5, and 1.0). These data cover a wide range of design 
variables used in practice. The results are shown in tables 1 and II and 
figures 3 to 9. In the discussion to follow, the hydrodynamic load and side 
leakage will be termed simply as load and flow, respectively. Using Patir 
and Cheng [8] shear stress factors, the frictional force on the journal sur- 
face was computed. The frictional forces at some eccentricity ratios were 
unrealistic. Hence they are not reported herein. 

In table I the performance characteristics of a finite journal bearing 
calculated from the present method of solution for smooth surface bearings 
(A « 6) are compared with a similar available solution [9]. For c • 0.8 
the agreement is not as good as for lower eccentricity ratios. 

A parametric study of the bearing characteristics is made in the 
following paragraphs. 

Effect of roughness parameter A . - The variation of loaa, attitude 
angle, and flow of a bearing having L/D = 1.0 and -y * 1 is shown as a 
function of roughness parameter in figures 3 to 5 for various eccentricity 
ratios and variance ratios. The load and flow increase but attitude angle 
decreases with increasing eccentricity ratio. This is a typical character- 
istic of a hydrodynamic oil journal bearing under steady-state conditions. 

When both journal and bearing have the same roughness structure (i.e., 
V^j B 0.5), load increases; flow decreases with increasing roughness 
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(i.e., for small A). The attitude angle remains more or less constant with 
roughness parameter for the same configuration. 

When the journal surface is rough and the bearing surface is smooth 
(i.e., > 1.0), is positive. This makes the last term on the 

right side of the Reynolds equation positive, thereby giving less pressure. 
Hence load decreases with increasing roughness. But for smaller values of 
A (A< 2), the load increases. This strange behavior of W^ for small 
A can be understood by looking carefully at the expression for if 

is plotted against A, it will be seen that increases with A 
to start with and then decreases as A is further increased and will ap- 
proach 0 when A is very large. At smaller values of A the number of 
contact points increases, and these points permit little or no flow in the 
direction of motion thereby increasing pressure and hence load. 

On the other hand, when the journal surface is smooth and the bearing 
is *"00911 (i.e., *•0), 4'^ is negative. We expect load behavior just 

opposite to that of .= 1.0. This, in fact, is shown in figure 3 for 
= 0 . Thus we get som? intermediate values of load when = 0.5 
if other parameters are kept constant. 

From figure 4 it can be seen that for small e (0.2), the attitude 
angle is not influenced much by A. For a higher value of c (0.8), the 
attitude angle for = 1.0 drops sharply at small A and becomes nega- 
tive when A = 1. 

The side leakage is plotted against roughness parameter for various e 
and V . in figure 5. When V„. = 0.5, the side leakage decreases with 

' j ' j 

increasing roughness. As the roughness is increased, these is a likelihood 
of more asperity contacts. This restricts the flow. Under this situation 
(V^j = 0.5), it may be noted that the load increases (fig. 3). When 
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Y « 1, Therefore when is 1.0 or 0, we expect a some- 

wfidt similar behavior of flow with load (flow is dependent on i ), The 
present results (fig. 5) show the similar trend. 

Effect of surface pattern parameter y . - The load and flow are shown 
(figs. 6 and 7) with respect to roughness parameter for various surface 
pattern parameters and eccentricity ratios. The is higher for longi- 
tudinally oriented surfaces (t > 1) and lower for transversely oriented sur- 
faces (y < 1) than for isotropic surfaces (y ■ 1). Again 6^ (AH, y) * 

^ (Ah, 1/y). As load is a function of both t and 6 and as both 

X ^ y 

(i and t are related to y in such a way, it is quite difficult to 

X y 

foresee the optimum values of load and flow for a particular y. From fig- 
ures 6 and 7 it is observed that for a small eccentricity ratio (e ■ O.J?), 
the maximum load and mitiimuin side leakage Q are given by y * b for 
L/D = 1.0 and . 0.5. For the same bearing configuration, when 
e = 0.6, the isotropic surface gives maximum load, but the flow is minimum 
for a surface having y = 6* 

Effect of length-to-diaineter ratio L/D . - In figures 8 and 9 the per- 
formance behavior for various L/U and c with * 0.5 and y = 1 is 
shown. As expected, a bearing with higher L/D and c gives higher load 
and lower flow q. In table II the load and side leakage for various L/0 
ratios and y with e ■ 0.2 is given. As mentioned earlier, p , ^ , 

X y 

and Y related in such a way that it is extremely difficult to predict 
which particular bearing configuration, so far as L/D and y are con- 
cerned, gives maximum load and minimum side leakage. In view of this, a 
large amount of data was generated and tabulated (table II). It can be seen 
from table II that, when both journal and bearing surfaces have Identical 
roughness structures, a bearing with L/D » 2.0 and y * 1 gives maximum 
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load, but a bearing with L/0 - 0.5 and y - 6 gives minimum side leakage. 
However, the flow in this case will not be minimum, since Q is non- 
dimensionalized with L/D. The minimum absolute side leakage q will, of 
course, be given by a bearing having L/0 ■ 0.5 and y - 6 under such an 
operating condition. 

CONCLUSIONS 

By using a flow simulation model developed by Patir and Cheng [5], 
tho steady-state performance behavior of a finite oil journal bearing using 
Reynolds boundary conditions was obtained. From this study and the results 
presented the following conclusions were drawn: 

1. When the roughness parameter A is 6, the results (smooth surface 
solution) obtained from the present method of solution are in reasonably 
good agreement with ttie similar available solution. 

2 . The hydrcdynaniic load increases and dimensionless side leakage de- 
creases with increasing surface roughness when both journal and bearing sur- 
faces have the same roughness structure. 

3. When the journal surface is very rough (A < 2) and the bearing 
surface is smooth, the load capacity and side leakage Q are higher for 
surfaces having the same roughness structure. 

4. when the bearing surface is very rough (A < 2) and the journal 
surface i? smooth, the load capacity and side leakage Q are lower for 
surfaces having the same roughness structure. 

5. The values of hydrodynamic load and side leakage Q for bearings 
having tne same roughness structure on both surfaces are in between the 
values for the bearings having rough journal surfaces and the values for 
bearings having rough bearing surfaces. 
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6. When both surfaces have the same roughness structure and when the 
bearing operates under lightly loaded conditions, a bearing with a longi- 
tudinally oriented surface pattern gives maximum load and side leakage, and 
a bearing with a transversely oriented surface pattern gives minimum load 
and side leakage. 

7. The contact load of a partially lubricated surface can sometimes be 
'comparable to hydrodynamic load. 


REFERENCES 

1. Greenwood, J. A.; and Williamson, J. B. P.: Contact of Nominally Flat 

Rough Surfaces. Proc. R. Soc. (London), vol. A295, no. 1442;, Dec. 

1966, pp. 300-319. 

2 . Whitehouse, D. J.; and Archard, J. F.: The Properties of Random Surfaces 

of Significance in Their Contact. Proc. k. Soc. (London), vol. A31b, 
no. 1524, Mar. 1970, pp. 97-121. 

3. Nayak, F. R.: Random Process Model of Rough Surfaces. J. Lubr. Technol., 

vol. 93, no. 3, July 1971, pp. 398-407. 

4. Greenwood, J. A.; and Tripp, J. H.: The Contact of Two Nominally Flat 

Rough Sui faces. Proc. Inst. Mecli. Eng. (London), vol. 185, 1970-71 , 
pp. 625-633. 

5. Patir, N.; and Cheng, H. S.: An Average Flow Model for Determining 

Effects of Th*'ee-Dimensional Roughness on Partial Hydrodynamic 
Lubrication. J. Lubr. Technol., vol. 100, no. 1, Jan. 1978, pp. 12-17. 

6. Tender, K.: Simulation of the Lubrication of Isotropically Rough Sur- 

faces. ASLE Trans., vol. 23, no. 3, July 1980, pp. 326-333. 


18 



7. Majumdar, B. C.; and Hamrock, B. J.: Effect of Surface Roughness on 
Hydrodynamic Bearings. NASA TH-B1711t 19B1. 

6. Patir, N.; and Cheng, H. S.: Application of Average Flow Nodel to 
Lubrication Between Rough Sliding Surfaces. J. Lubr. Technol., vol. 
101, no. 2, Aor. 1979, pp. 220-230. 

9. Pinkus, 0.; and Sternlicht, B.: Theory of Hydrodynamic Lubrication. 
McGraw Hill Book Co., Inc., 1961, p. 86. 


19 



lAHIf. 1. - COMPARISON OF HYDRODYNAMIC LOAD CAPACITY, 
AUITUDE ANGLE, AND SIDE LEAKAGE WITH AVAILABLE 
SMOOTH SURFACE BEARING lOLUTlON 
[L/D • 1.0, A ■ 6 , ■ O.S.] 


B 

W|| (load) 

(attitude angle) 

Q (side leakage) 

0.2 

0.08 (0.08)‘ 

73.6 (74) 

0.31 (0.32) 

.4 

.19 (.20) 

61.6 (62) 

.61 (.61) 

.6 

.45 (.44) 

48.6 (50) 

.93 (.94) 

.8 

1.46 (1.18) 

33.3 (36) 

1.29 (1.24) 


‘Values in parentheses are from Pinkus and Sternlicht 
[9] (table 4-1) p. 86. 









TABLL 11. - VARIATION OF HYDRODYNAMIC 
LOAD AND Slut ItAKAGE WITH VARIOUS 
LENGTH-T0-D1AM£TFR RATIOS AND 
SURFACE PATTERN PARAMETERS 
[c = 0.2, V^j . O.S.] 
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Figure 1. - Diagnm of an oil journal bearing with a 
greatly exaggerated clearance. 



Figure X - Variation of hydrodynamic load with roughness parameter 
fo( various eccentricity ratios and variance ratios. 
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Figure 4. • Vtrlition of attitude angle with roughness parameler tor 
various accentricity ratios and variance ratios. 



Figure 5. - Variation o( side leakage with roughness parameter tor 
various eccentriaty ratios and variance ratios. 




Figure 6. - virielion of hydrodynimic loed with roughness pire.neler 
lor viriru roughness pettern pererntters. 



Figure 7. - Verlitton of side leekege rith roughness pirenctcr lor 
Mrious roughness pittern piremeters. 



